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POISSONIAN STATISTICS IN THE EXTREMAL PROCESS OF 
BRANCHING BROWNIAN MOTION 1 

By Louis-Pierre Arguin 2 , Anton Bovier 3 and Nicola Kistler 4 

Universite de Montreal, Rheinische Friedrich-Wilhelms-Universitat Bonn 
and Rheinische Friedrich-Wilhelms-Universitat Bonn 

As a first step toward a characterization of the limiting extremal 
process of branching Brownian motion, we proved in a recent work 
[Comm. Pure Appl. Math. 64 (2011) 1647-1676] that, in the limit of 
large time t, extremal particles descend with overwhelming proba- 
bility from ancestors having split either within a distance of order 1 
from time 0, or within a distance of order 1 from time t. The result 
suggests that the extremal process of branching Brownian motion is 
a randomly shifted cluster point process. Here we put part of this 
picture on rigorous ground: we prove that the point process obtained 
by retaining only those extremal particles which are also maximal in- 
side the clusters converges in the limit of large t to a random shift of 
a Poisson point process with exponential density. The last section dis- 
cusses the Tidal Wave Conjecture by Lalley and Sellke [.Ann. Probab. 
15 (1987) 1052-1061] on the full limiting extremal process and its re- 
lation to the work of Chauvin and Rouault [Math. Nachr. 149 (1990) 
41-59] on branching Brownian motion with atypical displacement. 

1. Introduction. Branching Brownian motion (BBM) is a continuous- 
time Markov branching process which plays an important role in the theory 
of partial differential equations [4, 5, 23], in particle physics [24] and in the 
theory of disordered systems [6, 14]. It is also widely used in biology to 
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model the genealogies of evolving populations, the spread of advantageous 
genes, etc., [15, 19]. It is constructed as follows. 

Start with a standard Brownian motion (BM) (we will often refer to 
Brownian motions as "particles"), x(t), starting at 0. After an exponential 
random time, T, of mean 1, the BM splits into k independent BMs, inde- 
pendent of x and T, with probability p^, where XlfcLiPfc = 1j SfcLi ^Pk = 2 
and K = k{k — l)pk < oo. Each of these processes continues in the same 
way as first BM. Thus, after time t > 0, there will be n(t) BMs located at 
xi(t),... ,x n (t)(t), with n(t) being the random number of offspring generated 
up to that time [note that Kn(t) = e ]. 

An interesting link between BBM and partial differential equations was 
observed by McKean [23]: denote by 

(1.1) u(t,x) = F max xdt)<x 

-\<k<n{t) 

the law of the maximal displacement. Then, a renewal argument shows that 
u(t,x) solves the Kolmogorov-Petrovsky-Piscounov or Fisher [F-KPP] equa- 
tion, 



u t = -u xx + > v k u k - u, 



k=l 

1.2) 

ri, ifx>o, 

K ' J \0, if ar<0. 

The F-KPP equation admits traveling waves: there exists a unique solution 
satisfying 

(1.3) u(t,m(t) +x) — >u;(x) uniformly in x as t — >■ oo, 
with the centering term, the front of the wave, given by 

(1.4) m ( t ) = v ^_ ^ 

and a; (a;) the unique (up to translation) distribution function which solves 
the ordinary differential equation 

(1.5) \w xx + V2u x + uJ 2 — uj = 0. 

The leading order of the front has been established by Kolmogorov, Petro- 
vsky and Piscounov [20], whereas the logarithmic corrections have been ob- 
tained by Bramson [8], using the probabilistic representation given above. 

The limiting law of the maximal displacement has been studied intensely. 
Let 

n(t) 

(1.6) Z(t) = YjWM ~ x k(t)) exp -V2(V2t - x k (t)) 

k=l 
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denote the so-called derivative martingale. Lalley and Sellke [21] proved 
that Z(t) converges almost surely to a strictly positive random variable, Z, 
and established the integral representation 

(1.7) uj(x) = E[e- CZe '^ x ], 

with C > a constant. Thus the law of the maximum of BBM is a random 
shift of the Gumbel distribution. Moreover, it is known that 

(1.8) 1 -u(x) ~xe-^ x , x^+oo, 

where ~ means that the ratio of the terms converges to a positive constant; 
see, for example, Bramson [8] and Harris [17]. (There is emerging evidence 
that right-tails such as (1.8), manifestly different from those of the Gumbel, 
play an important role in a number of different fields, e.g., in models on 
spin glasses with logarithmic correlated potentials by Carpentier and Le 
Doussal [12], and Fyodorov and Bouchaud [16].) 

Contrary to the maximal displacement, very little is known on the full 
statistics of the extremal configurations (first-, second-, third-, etc., largest) 
in BBM. Such statistics are completely encoded in the extremal process, 
which is the random point measure associated to the collection of points 
shifted by the expectation of their maximum, that is, the point process 

n(t) 

(1.9) E(t) = J2^(t), Xi(t)=Xi(t)-m(t). 

i=l 

The key issue of interest is to characterize the limit of this process, as t "f oo. 
It can be shown that the limit of the point process exists using Bramson's 
analysis [7] on the convergence of solutions of the KPP equations with ap- 
propriate initial conditions [9, 18]. 

For given realization of the branching, the positions {xi(t)}i<. n u\ form 
a Gaussian process indexed by i G {1, . . . , n(t)} = T,t with correlations given 
by the genealogical distance 

(1.10) Qij(t) = sup{s < t:xi(s) = Xj(s)} 

(the time to first branching of the common ancestor). The information about 
the correlation structure of any subsets of particles in BBM is encoded in 
their genealogical distance. This applies, in particular, to the subset of ex- 
tremal particles, for which the following result was proved in [3]: with proba- 
bility tending to 1, branching can happen only at "very early times," smaller 
than r<2 with = 0(1) as t — > oo, or at times "very close" to the age of the 
system, namely greater than t — r g for r g = 0(1) as t — > oo. (The reason 
for this notation, in particular the use of the subscripts, will be explained 
below.) More precisely, denoting by St(-D) = {i S S ( :xj(i) G D} the set of 
particles in the subset m(t) + D, we have: 
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no branching 
Fig. 1. Genealogies of extremal particles. 

Theorem 1 [3]. For any compact Act, 
(1.11) lim sup F[3i,j€Z t (D):Q ij (t)€(r d ,t-r g )}=0. 

rd ' r 9^ 00 t>3max{r d ,r 9 } 

Figure 1 presents a graphical representation of the genealogies of extremal 
particles of BBM. 

Theorem 1 gives insight into the limiting extremal process of BBM. In fact, 
it suggests the following picture, which holds with overwhelming probability 
in the limit when first t'foo, and rj, r g — > oo after that. 

First, ancestries in the interval [0, cannot be ruled out: this regime 
generates the derivative martingale appearing in the work of Lalley and 
Sellke [21]. Moreover, since the ancestors of the extremal particles evolved 
independently for most of the time (namely in the interval \r^t — r g \), the 
extremal process must exhibit a structure similar to a Poisson process. Fi- 
nally, since ancestors over the period [t — r g ,t] also occur, it is natural to 
conjecture that small grapes of length at most r g = 0(1), that is, clusters of 
particles with very recent common ancestor, appear at the end of the time- 
interval. (According to this picture, the subscript in refers to derivative 
martingale, while that in r g stands for grape.) 

It is the purpose of this work to make part of this picture rigorous. In Sec- 
tion 2 we present our main result, which is proved in Section 3. In Section 4, 
we introduce a cluster point process, which we conjecture to correspond in 
the limit to the extremal process of BBM. We also discuss the cluster point 
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process in relation to the work of Chauvin and Rouault [13] on BBM con- 
ditioned to perform unusually large displacements, and in relation to the 
Tidal Wave Conjecture of Lalley and Sellke [21]. Detailed properties of this 
cluster point process will be the subject of a subsequent paper [2]. 

2. Main results. Despite the rather clear image described above, a frontal 
attack on the extremal process appears to be difficult. This is in particular 
due to the fact that one has to take into account the self-similarity of BBM 
which is first and foremost detectable in the small clusters, an issue which 
remains rather elusive (see Section 4 for more on this). On the other hand, 
the picture naturally suggests the existence of an underlying point process 
obtained from the extremal particles by a thinning procedure, which we 
describe next. 

Assume that the positions of particles at time t are ordered in decreasing 
order: 

(2.1) x 1 (t)>x 2 (t)>--->x n{t) (t). 

The inequalities will in fact be strict for almost all realizations of BBM for 
any deterministic time t. Define also 

( 2 - 2 ) Q{t) = {Qij W}i,j<n(t) = {*~ Qii(*)}i,j<n(t)- 

The pair (H(i), Q(t)) admits the following natural thinning. Since the matrix 
Q(t) is constructed from the branching of the BBM, the relation Qij(t) > q 
is transitive for any q > 0: 

(2.3) Qij(t)>q and Q jk (t)>q => Q ik (t) > q. 

In particular, for any q > 0, this relation defines an equivalence relation on 
the set {1, ... , n(t)}. The corresponding equivalence classes are just the par- 
ticles at time t that had a common ancestor at a time later than tq. We want 
to select a representative of each class, namely the maximal particle within 
each class, and then consider the point process of these representatives. For 
any q > 0, the q-thinning of the process (E(t),Q(t)), denoted by 3^(t), is 
defined recursively as follows: 

k = 1; 

(2-4) 

i k =mm{j > < <?, VZ < k- 1}; 

and 

(2.5) 3®(t) = (S$\t),k€ N) = (x ik (t), k G N), 

where it is understood that E)i^(t) = when an index i k in {1, . . . ,n(t)} 
satisfying min{j > z^-i -Qij < Q V/ < k — 1} can no longer be found. The 
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Fig. 2. Cluster- extrema. 



procedure selects the maximal position in each equivalence class defined from 
the relation Qij(t) >q. In addition, it is easily checked that the thinning 
map, 

(2.6) (H(t),Q(t))->Efo>(*), 

considered at the level of realizations, is a continuous function on the space of 
pairs (X, Q), where X is a sequence of ordered positions and Q is a symmet- 
ric matrix with entries in [0, 1], satisfying (2.3) (when this space is equipped 
with the product topology in each coordinate of X and Q). 

The thinning map can also be applied to t-dependent values of q. For ex- 
ample, take q = q(t) = 1 — r g /t, where r g is fixed t. In this case, the thinning 
effectively retains those particles which are extremal within the class defined 
by a "very recent" common ancestor, which we refer to as cluster- extrema. 
Figure 2 presents a graphical representation of the set of such particles. 

Our main result states that all such thinned processes converge to the 
same randomly shifted Poisson Point Process (PPP for short) with expo- 
nential density. 

Theorem 2. For any < q < 1, the processes z.( g \t) converge in law to 
the same limit, 2°. Also, 

(2.7) lim lim E^~ rB ^(t) = 2°. 

Moreover, conditionally on Z, the limit of the derivative martingale (1.6), 

(2.8) E° = PPP(C-Z-V2e-^ x dx), 
where C > is the constant appearing in (1.8). 
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The point process H has a fundamental connection with the limiting ex- 
tremal process of BBM. To see this, suppose for simplicity that the processes 
(E(t),Q(t)) induced by the law of BBM converge, as tfoo, to a process, 
(S, Q). (The laws of these processes are in fact tight because the law of S(t) 
is itself tight; see, e.g., Corollary 2.3 in [3], and that Qij(t) € [0, 1] for any 
Convergence would evidently follow from a complete characterization of the 
extremal process.) It follows from Theorem 1 that is either or 1. This 
suggests: 

(1) to define a cluster of particles as the maximal set of particles such that 
Qij = 1 for all i,j in the set; 

(2) to look at the process of the maxima of each cluster, denoted by, say, S°, 
defined as in (2.4), but where = min{j > i^-i : Q^j = V7 < k — 1}. 

We claim that H° is in fact the limit H° of E«)(t) in Theorem 2. Indeed, 
in view of the continuity of the thinning map (2.6), H^(i) converges to the 
g-thinned process, 3™, constructed from (H, Q) for all q. But, for any < 
q < 1, the g-thinned processes, H^, constructed from (S, Q) using (2.4) are 
equal trivially to S°, since is either or 1. The claim then follows from 
Theorem 2. The point process describing the particles at the frontier of BBM 
in the limit of large times is thus formed by two "types" of particles: those 
coming from the randomly shifted PPP with exponential density, the cluster- 
extrema; and the second type of particles, those forming the clusters. Clearly, 
particles in the same cluster always lie on the left of the corresponding 
Poissonian particles, by the very definition of the cluster-extrema. It remains 
an open question to characterize the law of the clusters (see Section 4 for 
some conjectures). 

We remark that, since r g = O(l) as t — > oo, the thinned process sj 1 r9// ^ 
is obtained from the extremal one by removing only a small number of 
particles, those which have genealogical distance smaller than t — r g from 
the maximum in their class. It is rather surprising at first sight (but not quite 
when seen under the light of Theorem 1) that such a point process converges, 
despite the high correlations among the branching Brownian particles, to 
a PPP with exponential density. 

Theorem 2 also provides insights into a result by Bovier and Kurkova [6] , 
who addressed the weak limit of the Gibbs measure of BBM, the random 
probability measure on £t attaching weights 



(2.9) 



exp(/fa fc (t)) 
UP) 




where f3 > is the inverse of temperature. To see this, let us first recall the 
following. 
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Consider the random set £ N) where the £'s are generated according 
to a PPP with density CZ\/2e~^ x dx on the real axis, C and Z as in 
Theorem 2. Construct then a new random set (pi,i G N) where /Oj = exp(/3£j). 
For /3 > v2, it is easily seen that A/"(p) = ■ Pj < oo almost surely, in which 
case the normalization pi = pi/N{p) is well defined, and the law of the 
normalized collection (pi,i £ N) is the Poisson-Dirichlet distribution with 
parameter m(/3) = \^2//3, which we shall denote by PD(m(/3)). 

In a somewhat indirect way (by means of the so-called Ghirlanda-Guerra 
identities, which avoid the need to first identify the limiting extremal pro- 
cess) Bovier and Kurkova proved that, in the low temperature regime (3 > 
y/2, the Gibbs measure Gp,t converges, in the limit of large times, to the 
PD(m(/3)); together with our Theorem 2, this naturally suggests that the 
Gibbs measure of BBM is concentrated, in fact, on the cluster- extrema. 

Finally, Theorem 2 sheds light on a property of the extremal process of 
BBM that was conjectured by Brunet and Derrida [10]. They suggested that 
the statistics of the leading particles are invariant under superposition in 
the sense that the extremal process of two independent branching Brownian 
motions has the same law, up to a random shift, as the extremal process 
of a single one. This property at the level of the entire process is likely to 
involve specific features of the laws of the individual clusters. On the other 
hand, at the level of the thinned process, it is a straightforward consequence 
of Theorem 2, since the law is Poisson with exponential density. 

Corollary 3. LetE(t) andr!{t) be the extremal processes (1.9) of two 
independent branching Brownian motions. Denote by Z and Z' the pointwise 
limit of their respective derivative martingale. Then, for any < q < 1, the 
law of the q-thinning o/H(i) conditionally on Z and Z' converges to 

(2.10) WP{C-{Z + Z')-V2e-^ x dx). 

In particular, the thinned process of H(i) + H'(t) has the same law in the 
limit as the thinned process H° of a single branching Brownian motion, up 
to a random shift. 

As mentioned before, Theorem 2 is a natural consequence of Theorem 1. 
The main ingredient is the following lemma, which allows to compare thin- 
ning processes on a set of large probability. We use the notation 

(2-11) (t ) | (^oe) ^ {S^ (t ) : (t ) > y} . 

Lemma 4. For any y£l and any e > 0, there exists ro = ro(y,e) such 
that for rd,r g > ro and t > 3max{r 5 , r^}, on a set of probability 1 — e, 

(2-12) H(")(t)| (j/i00) =H^ /4) | (y>00) , 
for any < q < 1 — r -j- . 
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Theorem 2 is then proved by a standard Poisson convergence argument 
which exploits the weak correlations between the cluster-extrema in classes 
of the ^-thinning. 

Proposition 5. With C > and Z the limiting derivative martingale, 
conditionally on Z, 

(2.13) lim lim E^*\t) = PPP(CZy/2e~ y/ * x dx). 

3. Proofs. 

Proof of Lemma 4. Theorem 1 describes the genealogies of particles 
which fall into compact sets around the level of the maximum but for the 
proof of Lemma 4 we need a slight extension in order to cover the case of 
sets which are only bounded from below; more precisely, we claim that for 

(3.1) lim sup F[3i,je-E t (y,oo):Q ij (t)e(r d ,t-r g )]=0. 

r d, r g-^°° i>3max{r d ,r 9 } 

To see this, we recall the following estimate proved by Bramson [8], Propo- 
sition 3: 

(3.2) P[ max sjfc(i) > y] < k(Y + l) 2 e~^ Y , 

-k<n(t) 

which is valid for t > 2,0 < Y < y/t and k > a numerical constant. The 
bound (3.2) implies in particular that 

(3.3) lim supP[(tS t (y, oo) > 0] = 0. 

y->oo t > 2 

For Y > y, using the splitting T, t (y, oo) = T, t (y,Y) U E t (Y,oo), we have the 
bound 

F[3i,jeZ t (y,oc):Q ij (t)£(r d ,t-r g )} 

<P[3i,j£ MV, Y) : Qij (t) € (r d , t-r g )}+ P[»S t (Y, oo) > 0] . 

The first term on the right-hand side vanishes, by Theorem 1, in the limit 
t — > oo first and r d , r g — > oo next, whereas the second term vanishes, by (3.3), 
in the limit t — > oo first and Y — > oo next: this proves (3.1). 
Let us denote by A t ^ d ,r g (y, e) the event 

(3.5) {3i,j e E t (y,oo):Qij(t) G [r d ,t-r g ]}. 

By (3.1), there exists ro = ro(y, e) such that, for r d , r g > ro, P[A^ ] > 1 — e. 
By definition, on the event A^ r r , the following equivalence hold's for any 
??< q <l- r f: 



(3.6) Qij(t)<q ^ Qij(t)< 



t 

t ' 

The assertion of the lemma is now a direct consequence of the definition of 
the thinning E^(t) in (2.4). □ 
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To prove Proposition 5, we need some control on the derivative martingale. 
Lemma 6. Let 

(3.7) Z^{t) = {v / 2t-x fc (t)} 2 exp[-2v / 2{v / 2t-x fc (t)}]. 

k<n(t) 

For any given e > 0, 

(3.8) lim F[Z {2 \t) >e] =0. 

t— >oo 

Proof. First, by Bramson's estimate [8], we may find Y = Y(e) large 
enough, s.t. 



(3.9) 



maxx fe (t) - m(t) >y] < (1 + Y) 2 e~^ Y < e/2. 

k 



Using this bound, and the Markov inequality, we get 

F[Z^(t)>e] 
(3.10) t 

< —E[{V2t - x^lV 2 ^^*-^;^*) < m(t) + Y] + e/2. 



The first term on the right-hand side is bounded from above by 



e 

e 



t r°° 



(3/(2 v / 2))logt-Y 



1 

< - 

s 



xV 2 ^exp 



(V2t-x) 2 \ dx 



2t 



(3.11) < 



(3/(2v / 2))logt-y 

exp -v^(3/(2v/2) logt - Y) 



dx 
/2vft 



(3/(2 v / 2))logt-y 



dx 
J2wt 



< 



, . £-3/2 poo 



x 2 e~ x2 / 2t > 



dx 
/2tt7 



< 



p-t 



-3/2 



-t — > as i — > oo. 



This proves the lemma. □ 

Proof of Proposition 5 . We will show the convergence of the Laplace 
functionals. For 4> : R — > M + measurable with compact support, we claim that 



(3.12) 



lim lim E 

r c i~- >co t—>-oo 



exp 



^E^^t^dx) 



exp-CZ J (1 - e'^y^e'^ dx 



from which the proposition would evidently follow. 
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We will prove (3.12) for simple step functions, that is, of the form <p(x) = 
^2i=i a i^-Ai for a« > 0,i = 1, . . . ,N and Ai = [Aj, Aj],i = 1 . . . N disjoint com- 
pact subsets. The extension to the general case of measurable follows by 
a standard monotone class argument. 

We will make use of the splitting 

(3.13) m(t) = V2r d + m(t - r d ) + R t 
for some Rt = o(l) as 1 1 oo. 

We also introduce, for j < n(r d ), independent BBMs {x^\t — r d )} k<n .( t -r d ) > 
and use the abbreviation 

(3.14) Mjit-ra) 



max xi^(t — r d ) — m{t — r d ). 

k<rij(t-r d ) 



Conditionally on everything that happened up to time r d , the following 
equality holds in law due to the Markov property and the definition of the 
extrema in the (r^/^-thinning class: 

(3.15) 3('</*>(t) ( = } { Xj (r d ) - V2r d + M 3 {t - r d ) + R t } j=1 . Mrd y 

Since the Mf's are i.i.d., with ~^M(t—r d ) standing for expectation with respect 
to M(t-r d ), 



(3.16) 



As t — > oo, the variable M(t — r d ) converges weakly to M with law u by (1.3). 
Hence 



E 


exp — 


j cj){x)^ rdlt) {t){dx) 








M r d) 






= E 




(r d )-V2r d +M(t-r d )+R t ) 






■ i=i 





(3.17) 



lim E 

t— >oo 



exp 



:{r d /t) 



(t)(dx) 



E 



(3.18) 



Define Uj(r d ) = \f2r d — Xj(r d ). We write 

EA/[e~^ ( " w(rd)+M) ] = l-E A/ [l-e-^ ( -^ (rd)+A/) ] 
= :l-F(- yj (r d )), 
and 

(3.19) (3.17)=Eexp| log[l - F(- yj (r d ))] 

j<n(r d ) 
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Observe that 
(3.20) 



min Mj(rd) too 

j<n{r d ) 



a.s. 



as r d t oo. This implies that 



(3.21) 



max F(-y,j(r d )) |0. 

j'<n(r d ) 



Using that — rc — x 2 < log(l — x) < — x for < x < 1/2, for large enough, 
we obtain (up to a vanishing error) upper and lower bounds of the form 



E 



expj- ]T F(-y,(r d ))j 



(3.22) 



j'<n(r d ) 

> (3.19) 

> E 



expj- ^ F(- yj (r d )) - F(- yj (r d )f 

L ^ j<n(r d ) 

Since (ft is chosen to be a simple step function, 

(3.23) F (_ y .( rd )) = ^(i_ e -^) / da;. 



i=l 



Hence we can make use of the asymptotics (1.8) to obtain 



N 



n-Vj(r d )) ~ J> - e-^M^ + y.MK 



--^(Ai+i/jCrd)) 



i=l 



(3.24) 

-(^ + yj (r d ))e-^ + %^»}, 

with ~ meaning that the ratio of the left- and right-hand sides converges 
to 1, in the limit r d — > oo, P-a.s. We regroup the terms on the right-hand 
side to get 

F(- yj (r d ))^C yj (r d )e-^^ 

(3.25) 

N 

-a i){e -V2A i _ e -V2A i}+n{y .( rd))j 



i=l 

with 1Z containing all the remaining terms; clearly, 
(3.26) \n( yj (r d ))\<p-e-^"<r*\ 

where p depends on the Oj and Ai, but not on Uj(r d ). By the convergence 
of the derivative martingale as r d t oo [cf. (1.6)], and the fact that, in the 
same limit, 

j<n(r d ) 



(3.27) 
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/v 



(3.28) 



lim V F(x J (r d )-V2r d ) = CZT(l-e- a *){e 

r d t°o * — ' * — ' 
j<n{r d ) 1=1 



2-1, 



CZ / (1 



-<j>{x) 



)V2e 



dx, 



P-almost surely. This yields the correct asymptotics for the upper bound 
in 3.22. 

The lower bound in (3.22) involves exactly the same term as the left-hand 
side of (3.28), and the additional term 

(3.29) Yl F(xj(r d )-V2r d ) 2 . 

j<n{r d ) 

It is straightforward to see that (3.29) converges to zero, as r d too. In fact, 
by the same argument as in (3.23)^(3.26), one sees that 



r d too. 



(3.30) |(3.29)| =0\ 
With the notation of Lemma 6, 

(3.31) Y yAr d ?z- 2 ^ {Td) = Z [2 \r d ), 

j<n{r d ) 

and this converges to zero in probability, by Lemma 6. Hence, in the limit 
of large r d , the lower and upper bounds in (3.22) coincide, which concludes 
the proof of the proposition. □ 



Proof of Theorem 2. Let 
support. We need to show that 



be measurable, with compact 



(3.32) 



lim E 

t— >oo 



exp ■ 



:E 



exp-CZ J (l-e-^v^e-^dx 



for any ^ < q < 1 



.-. This is straightforward in view of Lemma 4 and 
Proposition 5 by taking y smaller than the minimum of the support of eft 
and e arbitrarily small. □ 



4. Open problems. 

4.1. On the extremal process of BBM. We consider the following cluster 
point process. Let (Q', J 7 ' , P), C > be a probability space, and Z : Cl f — > M + 
with distribution as in Theorem 2. (Expectation w.r.t. P will be denoted 
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by E.) Conditionally on a realization of Z, let (rjf,i G N) be the position of 
particles generated according to a Poisson point process with density 

(4.1) CZ(-xe-^ x )dx 

on the negative axis. For each i € N, consider independent branching Brow- 
nian motions with drift —y/2, that is, {ail (r) — \/2r;k < nj(r)}, issued on 
(fi', J-' ,P). ("Time" is denoted here by r.) 
Remark that for given i 6 N, 

(4.2) max x¥\r) — \Ph — > — oo, 

fc<Tij(r) 

P-almost surely. The branching Brownian motions with drift are then su- 
perimposed on the Poissonian points, that is, the cluster point process is 
given by 

(4.3) U r = {7r itk (r);ien,k = l...n i (r)}, Tr i>k (r) = m + xf{r) - V2r. 

The existence of the large time limit of Il r is not straightforward. Due 
to (4.2), only those Poissonian points whose attached branching Brown- 
ian motion performs an unusually large displacement can contribute to the 
limiting object. It is thus not clear that one finds any Poissonian points 
at all which, together with their cluster of particles, achieve this feat. The 
fundamental observation here is that, in virtue of (4.1), the density of the 
Poissonian points on the negative axis grows (slightly faster than) exponen- 
tially when x — > — oo. Together with the work of Chauvin and Rouault [13] 
on branching Brownian motions conditioned to perform unusually large dis- 
placements, this observation can be exploited to rigorously establish the 
existence of the point process n r in the limit of large times, as well as some 
of its statistical properties. We will report on this in a subsequent paper [2]. 

Here, we only put forward the following conjecture, which appears rather 
natural in the light of Theorem 1 and the results on the paths of extremal 
particles in BBM established in [3]: 

Conjecture 7. In the limit of large times, the distribution of the ex- 
tremal process of BBM, H(i) and that of II r coincide, that is, 

(4.4) lim 3(t) = lim II r . 

t— >oo r—¥oo 

In particular, with </>:R — > M + a measurable function with compact support, 



lim E 

t— >oo 

k<n(t) 

(4.5) 

fO 



exp (~ X] 4>( x k(t) -m(t))\ 

^ k<n(t) ' 
rO 

exp-CZ (l-e-^(*)){- xe -^ x }dx 

J — oo 



lim E 
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where 

(4.6) e"^ = E 



expf- 0(x + Xfc(r) - V2r) ) 



k<n(r) 



We remark that densities of the form — xexp(— \/2x) dx on the negative 
axis have been conjectured to play an important role in the recent work by 
Brunet and Derrida [10], where the average size of the gaps between the reth- 
and (re + l)th-leading particle at the edge of BBM is numerically shown to 
behave as 

(4.7) !__!_ + ..., 

n n log n 

(which is indeed "close" to the average size of the gaps in a PPP with density 
— xe~^ x dx on the negative axis). 

4.2. On a conjecture by Lalley and Sellke. Conjecture 7 is similar but 
fundamentally different from the Tidal Wave Conjecture formulated by Lal- 
ley and Sellke [21]. Lalley and Sellke suggested that the Poisson point process 
entering into the construction of IL. should have density CZe~^ x dx con- 
ditionally on a realization of Z where C is some constant. However, this 
cannot be correct. We will show that such a point process does not exist in 
the limit r — > oo: the density of the Poissonian component cannot compen- 
sate (4.2) and all the particles are bound to drift off to — oo. To formulate 
this precisely, consider the point process 

(4.8) ft, = (fji + xf (r) - \/2r; i G N, k = 1, . . . , n^r)), 

where the r/'s are points of a PPP with density CZe~^ x dx, and the x^'s 
independent BBMs. 

Proposition 8. For given j/gR, 

(4.9) lim P[U r [y, oo) >1\Z] = 0. 

1 — >oo 

In order to prove Proposition 8, we make use of the following bound 
established by Bramson: 

PROPOSITION 9 ([7], Proposition 8.2). Let yo<0 (strictly). There exists 
r = r o(yo) such that for r > ro, x > m(r) + 1 and z = x — rre(r), 

r i f0 p -(x-y) 2 /2r 

(4.10) P maxx k (r)>x < p ■ e r / (1 - e- 2{y - yo)z/r ) dy, 



lk<n(r) 



where p> is a numerical constant. 

Using this with yo := —1, we obtain the following corollary. (Here and be- 
low, p > denotes a numerical constant, not necessarily the same at different 
occurrences.) 
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Corollary 10. For X > 1, and r>r a = r D (— 1), 
P 



(4.11) 



max Xk(r) — m(r) > X 

k<n(t) W W - . 



<p.X.e W (-V2X-^ + ^X l0gr 



2r 2\/2 

Proof. According to Proposition 9, for X > 1, 
P 



(4.12) 



max Xfc(r) — m(r) > X 

■k<n(t) 



<p-e r 



e -(X+m(r)-j/) 2 /2r 

.1 V27rr 



(1-e 



-2(2,+1)X/tm 



dy. 



Since y + 1 > we have that 1 - e -Hy+^)^/r < 2(y + l)X/r. Using this, the 
right-hand side of (4.12) is at most 



(4.13) 



Xe r 



e _(x+m(r)-j/) 2 /2r 

(y + 1) tf= dy. 



Expanding the square in the Gaussian density, (4.13) is at most 



X 2 3Xloer 



p-X-exp -V2X + 

2r 2y/2r 



(4.14) 



Xy/r+V2y+y{3/(2V2))(logr)/r-y 2 /2r 



<i 



<p.X.expf-V2X-^ + ^^ 
V 2r 2\/2 r 

settling the proof of the corollary. □ 

PROOF OF Proposition 8. In view of (4.2), it is plain that for any 
finite set I C N 



(4.15) 



max 



|?7j+ max [x$ (r) — y/2r] \ — 



-OO, 



P-almost surely. But the number of Poissonian points (rji;i £ N) in the in- 
terval [0,oo) is finite, P-almost surely: this follows from the fact that the 
density CZe~^ x dx is integrable on x £ [0,oo). Hence, Proposition 8 will 
follow as soon as we prove that 



(4.16) P\3 ien :f]i + max {x$(r) - V2r} > y and fji G (-oo,0)|Z 



k<rii(r) 



rfoo 



0. 
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By the Markov inequality, and using that the BBMs superimposed on the 
Poissonian points are identically distributed, (4.16) is at most 



(4.17) 



P 



max {xk(r) — V2r} > y — x 

k<n(r) 



CZe'^dx. 



We rewrite this in terms of M(r) = maxj< n ( r ){xt(r) — m(r)}: 



(4.17) 



(4.18) 



P 



M(r) >y — x + 



2V2 



logr 



CZe-^dx 



{CZe 



-V2y\ 



P[M(r)>X]e^ x dX, 



r 3/2 Jy+3/(2V2)\ogr' 

the last step by change of variable y — x + logr — > X. 

Let us abbreviate p = CZe~^ 2y . (Note that y and Z are fixed.) For 
large enough, 

3 

(4.19) y + :r ^=logr> 1, 



hence we may use (4.14) to get that (4.18) is at most 

p r xex P (^x^y* 2 /(2r )dx 

'l/+(3/(2V^))logr \2V2 r J 

9 (logr) 2 ' 



(4.20) 



r 3/2 



P 

,.3/2 



cxp 
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=l+o(l), rfoo 
3 



by change of variable X — ^= log r — > Y . 
It thus remains to control the term 



(4.21) 



r 3 / 2 A 2^ r 3 / 2 L 



As for the first term on the right-hand side of (4.21): 



(4.22) 



1 POO 1 POO 

-4tz Ye~ Y ^ 2r dY = ^ xe~ x2 / 2 dx^0 as rfoo. 
r 3 / 2 J y J y/V¥ 



As for the second term on the right-hand side of (4.21): 
3 logr ro ° 



(4.23) 



2\/2 r 3 / 2 

3 logr 
~ 2V2 r 



e~ Y ^ 2r dY 



/ e~ x ' 2 dx -> as r 1 00. 



This proves (4.16), settling Proposition 8. □ 
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Remark 11. The above computations also suggest that a point process 
which is obtained by superimposing independent BBMs with drift —y/2 on 
a PPP with a certain density exists in the limit of large times if and only 
if such density is, up to a (possibly random) constant, — xexp (-y/2x) dx on 
the negative axis. 

In fact, a closer look at the above considerations reveals that the left-hand 
side of (4.22) is the leading order of the expected number of points (of the 
superimposed point process) which fall into the subset [y,oo). Choosing the 
density of the Poissonian component as in Conjecture 7, (4.22) would then 
read r _3//2 Y 2 e~ Y l 2r dY, which indeed remains of order 1 in the limit 
r — > oo. 

Note added in revision. There has been considerable activity concerning 
the extremal process of BBM after this paper was submitted for publica- 
tion. Brunet and Derrida have shown in [11] that all statistical properties of 
the rightmost points can be extracted from the traveling wave solutions of 
the Fisher-KPP equation. The validity of Conjecture 7 has been settled in 
a paper of ours [2] , where it is proved that the extremal process of branching 
Brownian motion weakly converges in the limit of large times to a Poisson 
cluster process; shortly after that, Aidekon et al. [1] recovered the same 
results by means of "spine techniques." The Poissonian structure of the ex- 
tremal process can also be proved using the property of superpos ability as 
observed by Maillard [22]. This property of the process was conjectured by 
Brunet and Derrida in [11] and proved in [2]. 

Acknowledgments. The authors thank Eric Brunet and Zakhar Kabluchko 
for interesting discussions on the existence of the extremal process of branch- 
ing Brownian motion. 
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